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I. Introduction 

Many real-world engineering systems are too complex to be defined in precise terms, imprecision is often involved in any 
engineering design process. Fuzzy systems have an essential role in this fuzzy modeling, which can formulate uncertainty in 
actual environment. In many linear systems, some of the system parameters are vague or imprecise, and fuzzy mathematics is 
a better tool than crisp mathematics for modeling these problems, and hence solving a fuzzy linear system [11] or a fuzzy 
differential equation is becoming more important. The concept of fuzzy numbers and arithmetic operations with these 
numbers were first introduced and investigated by Zadeh [22], A different approach to fuzzy numbers and the structure of 
fuzzy number spaces was given by Puri and Ralescu [16], Goetschell and Voxman [12] and Wu and Ma Ming [20], 

Since Friedman et al. proposed a general model for solving an FIX FI fuzzy linear systems whose coefficients matrix is crisp 
and the right-hand side is an arbitrary fuzzy numbers vector by an embedding approach, many works have been done about 
how to deal with some advanced fuzzy linear systems such as dual fuzzy linear systems (DFLS), general fuzzy linear systems 
(GFLS), full fuzzy linear systems (FFLS), dual full fuzzy linear systems (DFFLS) and general dual fuzzy linear systems 
(GDFLS). These works were performed mainly by Allahviranloo et al. [2, 3,4,5], Abbasbandy et al. [1,7], Zheng et al. [1,9] 
and Dehgham et al. [10] and so on. In traditional fuzzy linear systems, the uncertain elements were usually denoted by the 
parametric form of fuzzy numbers. Based on arithmetic operations of the fuzzy number, the fuzzy linear systems could be 
extended into crisp function linear systems. Therefore the solutions of the fuzzy linear systems can be obtained by solving 
the model by means of ordinary analytical and numerical methods. 

However, very few researchers have developed methods to solve fuzzy complex system of linear equations. The fuzzy 
complex numbers was introduced by J.J. Buckley in 1989 [9]. The ft X ft fuzzy complex linear systems have been studied 
by M.A. Jahanigh [15]. Solution of fuzzy complex linear system of linear equations was described by Rahgooy et al. [18] and 
applied to circuit analysis problem. In 2014, Behera and Chakraverty [8] discussed the fuzzy complex system of linear 
equations by adding and subtracting the left and right bounds of the fuzzy complex unknowns and the right-hand side is 
fuzzy complex vector. 

In this paper complex fuzzy linear system is investigated. A numerical procedure for calculating the fuzzy solution is 
designed. Finally, some examples are given to illustrate our method. The structure of this paper is organized as follows: 

In Section 2, we recall the triangle fuzzy number and present the concept of the complex fuzzy linear systems. The 
computing model to the complex fuzzy linear systems is proposed in detail and the fuzzy approximate solution of the 
complex fuzzy linear systems is obtained by solving the crisp systems of linear equations using QR- decomposition method. 
Some examples are given to illustrate our method in Section 4 and the conclusion is drawn in Section 5. 

II. Preliminaries 

There are several definitions for the concept of fuzzy numbers. 

2.1 The fuzzy number 
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Definition 2.1. A fuzzy number is a fuzzy set like U : R — > I = [0, l] which satisfies: 

(1) U is upper semi-continuous, 

(2) U is fuzzy convex, i.e., u(Ax + (l — /l) y) > min{tt(x), u(y ) } for all X,y R, A <E [0, l] , 

(3) U is normal, i.e., there exists X 0 e R such that u(x (t ) = 1 , 

(4) Slip pu = [x e R | u(x) > 0} is the support of the ll , and its closure clfsupp U) is compact. 

Let E be the set of all fuzzy numbers on R . 

Definition 2.2. A fuzzy number u in parametric form is a pair|//,//j of functions u ( /') , U ( r) , 0 < r < 1 ;which satisfies 
the requirements: 

(1) ll ( /' ) i s a bounded monotonic increasing left continuous function, 

(2) U ( r ) is a bounded monotonic decreasing left continuous function, 

(3) w(r)<w(r),0<r <1. 

A crisp number x is simply represented by (w ( r) , U ( r )) = (x, x), 0 < r < 1 . By appropriate definitions the fuzzy number 

space {u (r) ,u ( /' ) ( becomes a convex cone E ] which could be embedded isomorphically and isometrically into a Banach 
space. 

Definition 2.3 Let x = fx(r),x(r)V y = ^ y (r), y (r)j,0 < r < 1, and k e R . Then 


(1) 

t = T iff x(r) = y(r),x( 

r ) = y( r ). 

(2) 

x + y = (A(r) + y(r),x(r) + } ; (r)), 

(3) 

1 

1 

II 

1 

K 

-z( r ))’ 

(4) 

[kx,kx\,k > 0 
kx = < , 

( kx,kxj,k < 0 



Definition 2.4 An arbitrary fuzzy complex number may be represented as X = p + iq , where 
p = and q = r) ,^(r) j for all 0<r<l are two real fuzzy number. Hence one may have 

X =(p(r),p(r)) + i(q(r),q(r))■ 

I)efimtion 2.5 For any two arbitrary complex fuzzy numbers X = p + iq, y = U + iv where p,q,U,V are fuzzy numbers, 
their arithmetic is as follows: 

(1) z l + z 2 =(x l + x 2 ) + i(y l + y 2 ); 

(2) kx = kp + ikq, k e R\ 
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(3) Z t x z 2 - {(x 2 xx 2 )-(y, X y 2 )} + i {(x, X y 2 ) + (j 1 xx 2 )}. 

2.2 Complex fuzzy linear systems 

Definition 2.6. The linear system equation 

GA +c l2 Z 2 +... + c ln z n =w, 
c 21 z l +c 21 z 2 +- + c 2 n z„ =w 2 

< 

C nlZl+C n 2 Z 2 +... + C nn Z n =W n 


where C-, 1 < i, j < n are complex numbers and VtJ , 1 < i, j < n are complex fuzzy numbers, is called a complex fuzzy 
linear system(CFLS). 

Using matrix notation, we have 

Cz = W (2.2) 

A complex fuzzy numbers vector 
■c — (~1 ’ ^2’ Z n ) 


is called a fuzzy solution of the complex fuzzy linear system (2.1) if z satisfies (2.2). 

III. Solving complex fuzzy linear system of equation 

3.1 Equivalent fuzzy linear system 

Theorem3.1. The n x n complex fuzzy linear system (2.1) is equivalent to a 2/7 X 2/7 order fuzzy linear system 

Gx = b 


where 



( A - B ' 


M 



G = 

k b A , 

,x = 


,Z/ = 

v v y 


(3.1) 


(3.2) 


Proof. We denote C = A + iB , A, B e R' and W = U + iv where it , V are fuzzy number vectors. We also suppose the 
unknown vector z = p + iq where p, q are two unknown fuzzy number vectors. 

Since Cz = W , we have 

(A + iB)(p + iq ) = u + iv 


(. Ap - Bq ) + i{Aq + Bp) = u+ iv . 

Comparing with the coefficient of i , we have 

f Ap -Bq = u 
[. Aq + Bp = v 


'A -B" 

V 



A y 



Ay 
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It admits a 2 n order fuzzy linear system. 

We express it in matrix form as follow: 

Gx = b 

3.2 The model 

When z, W of fuzzy linear equation (2.3) are denoted by the parametric form i.e. 

w = [w x , w 2 , ..., w n f , Wj = Uj + ivj = [u j (r),u j (r)\ + i[Vj (r ) , vj (r )], 

Z = [z v Z 2 ,...,Z n ] T ,Zj= Pj+iqj = [p.(r),p j (r)] + i[q,(r),q j (r)], j = l,2,...,n,0 < r < 1. 

Some results for solving fuzzy linear equation (3.2) are obtained by the following analysis. 

Theorem 3.2. The fuzzy linear equation (3.2) can be extended to a crisp function linear system as follow 

SX(r) = Y(r) 

Where 


5 = 

r G + -G A 

,X(r) = 

r x(r)^ 

,F(r) = 

r b(r ) ^ 


v-G 


{~ x ( r )y 


y-H r )y 


(3.3) 


(3.4) 


in which the elements gt of matrix G + and gy of matrix G + are determined by the following way: if 
gy ^ 0, g+ = g 0 else gy =0,1 <i,j< n; if gy < 0, g~ = g y else g~ = 0, 1 < i, j < n. 

Proof. Let b = [Z?(r),Z?(r)],0 < r < 1 and x = [x(r),x(r)]. Suppose G — G +G in which the elements gy of 
matrix G + and gy of matrix G + are determined by the following way: if gy > 0, gy = gy else gt = 0, l</,;'<n;if 
gy < 0, gy = gy else gy = 0, 1 < /, j < Tl. The fuzzy linear equation Gx = b can be expressed as 

(G + +G}[x(r),x(r)] = [b(r),b(r)]. 

\{kxj{r),kxj{r)^,k> 0 


Since 


We have 


kXj=\ 


Gx = < 


( kxj (r^,kxy (r)),£ < 0 

(Gx(r),Gx(r)j,G > 0 
| Gx(r),Gx(r)j,G < 0. 


So the equation (3.4) be rewritten as 

G + [x(r),x(r)] + G“[x(r),x(r)] = [G + x(r),G + x(r)] + [G”x(r),G”x(r)], 
[G + x(r) + G“x(r),G + x(r) + G”x(V)] = [b(r^,b (r)]. 


Thus we have 
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jG + x(r) + G x(r) = Z?(r) 
jG + x(r) + G~x(r) = Z?(r) 

f G + x(r)-G”(-x(r)) =b(r) 

or < _ _ 

[G + (-x(r))-G~x(r) = —b(r) 

Expressing in matrix form , we have 


' G + -G " 

' *( r ) ' 


' b(r ) N 

v-G G\ 

~ x ( r \ 


-H r )y 


Remark 3.1 If Cz = W , then z + Z and z — Z are the solutions of complex fuzzy linear system, i.e., 

C(z + z) = (w+w\ and C(z~z) = (w-VtY 

Remark 3.2 If Cz = W, then Z~ Z is the solution of complex fuzzy linear system, i.e., c(z~z\ = ^ tv — H’j and 

|C|fz-zj = (w-wj. 

In order to solve the fuzzy system of linear equation (3.2), we need to consider the systems of linear equations (3.4). It seems 
that we have obtained the minimal solution of the fuzzy linear system (3.2) as 


Z(r) = 5 t F(r) 


i.e., 


^x(r)^ 


' G + -G" Y 

f -_ (r) ] 

A r )j 


v-G G + , 

-*(>), 


(3.5) 

(3.6) 


where S * is the Moore-Penrose generalized inverse of matrix S . 


Definition 3.1. Let X (r) = (xz (r),— Xj (r)Vl < j <2n denotes the minimal solution of (2.5). The fuzzy number 
vector Z = [p (r), p j(r)\ + i[q (r),q y ( r )] , 1 < j <n defined by 


Pj(r) = mm{Xj(r),Xj(r),Xj(i),Xj(i)}, 

Pj(r) = max{x . (r) , xj (r) , x . (l) , Xj (l)} , (3.7) 

j = l,2,...,n,0 < r < 1, 


q. ( r ) = minlx^. (r) , xj (r) , x . (l) , xj (1)} , 

q, (r) = max{x . (r),xj ( r),x } (l),xj (l)}, (3.8) 

j = n + \,n + 2,. ..,2/t,0 < r < 1. 

is called the fuzzy minimal solution of the fuzzy linear systems (2.3). If p (r ) , q . ( r ) , 1 < j < 2 n are all fuzzy numbers 

then z = { [P (r),pj (r)] + i[q ( r),q j (r)],l < j < n) is called a strong complex fuzzy minimal solution of the 
complex fuzzy linear systems (2.1) . Otherwise, Z is called weak complex fuzzy minimal solution. 
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3.3 QR-decomposition method 

Theorem 3.3 If eq.(3.4) S e in x n , then there exist an orthogonal Q e mxn and an upper triangular R e in x n so that 
S = QR. 

Proof. Suppose n = 1 and that Q is Householder matrix so that if R = Q T S then R (2 : in ) = 0 .It follow that S = QR 
is a QR factorization of S . For general n we partition S , 

5= [5, |v] 

Where V = S /i) . 

By induction, there exists an orthogonal Q 1 e R " so that *i = Qls l is upper triangular. Set W = Q' v and let 
H' ( n : in ) = Q 2 R n be the QR factorization of W ( n : in ) .If 


Q = Q , 



0 


a 


Then 

is a QR factorization of S . 


a = q[r 1 


w(l:n-l) 


] 


Lemma 3.1 If A is an in x n matrix with full column rank, then A can be factored as A = QR , where Q is an in x n 

matrix whose column vectors form an orthonormal basis for the column space of A and R is a nxn invertible upper 
triangular matrix. 

Theorem 3.4. If Eq.(2.5) S is a matrix with full column rank , and if S = QR is a QR — decomposition of A , then 
the system for SX — Y can be expressed as 

RX=Q t Y (3.11) 

and the solution can be expressed as 

X=R x Q t Y (3.12) 

Proof. According to Lemma 3.1., we have 

S = QR 


Where Q is 4m X 4/7 orthogonal matrix , R is a 4 n X 4/Z invertible upper triangular matrix . The fact that Q has 
orthonramal column implies that Q T Q = I , so multiplying both side of S — QR by on the left side 

R = Q t S (3.13) 


Therefore, the system for SX = Y can be expressed as 

RX = Q t Y . 

In order to obtain the solution, due to R is invertible upper triangular matrix, so 

X=R l Q T Y . 
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3.4 A sufficient condition for strong fuzzy solution 

The key point to make the solution matrix being a strong fuzzy solution is that S Y (r'j is fuzzy matrix , i.e., each element 
in which is a triangular fuzzy number. By the analysis, it is equivalent to the condition S' > 0 . 

Theorem 3.5. If 

((g + - G )' + (G + + G ) r ) > 0, ((g + - G ) ’ - (G + + G ) j ) > 0, 

the complex fuzzy linear equation(2.1) has a strong complex fuzzy minimal solution as follows: 

Z = [p j (r),p j (r)] + i[q,(r),q j (r)\,0<r<l (3.10) 


x(r) = 


f P{r)^ 


Where 


: ( r ) = 


v^( r ) j 

r p( r f 

v^( r )y 


= Eb(r') — Fb(r'} 

= -Fb(r^) + Ebi^r^ 


£ = I((G + -G ) f +(G + + G ) f ) 
F = i((G + -G’) t -(G + +G') t ) 


Proof. The proof is straight forward 

IV. Numerical Example 

In this section, we will demonstrate the efficiency and superiority of the proposed method using numerical examples. 
Example 4.1 Consider the following CFSLE: 


Z,~Z 2 =(r,2-r) + i(l + r,3-r) 

< ' 2 _ (4.1) 

Zj + 3 z 2 = (4 + r, 1 - 2r) + i(r - 4, -1 - 2 r) 


Let z l = p x + iq x = (p l ,q l ) + i(q l ,q l ), z 2 — P 2 + iq 2 = (p 2 ,q 2 ) + i(q 2 ,q 2 ). The extend maxtrix is 


r i 

o 

o 

1 

o 

o 

o 

CG 

i 

3 

o 

o 

o 

o 

o 

o 

o 

1 

1 

o 

o 

o 

o 

o 

o 

o 

1 

3 

o 

o 

o 

o 

o 

o 

o 

o 

1 

o 

o 

1 

o 

o 

o 

o 

1 

3 

o 

o 

o 

o 

o 

o 

o 

1 

1 

o 

vO 

o 

o 

o 

o 

o 

1 

3y 


And Rank(S)=8,so we obtain 
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.7071 

0.6396 

-0.2023 

-0.2236 

0 

0 

0 

0 

.7071 

-0.6396 

0.2023 

0.2236 

0 

0 

0 

0 

0 

-0.4264 

-0.6068 

-0.6708 

0 

0 

0 

0 

0 

0 

-0.7416 

0.6708 

0 

0 

0 

0 

0 

0 

0 

0 

-0.7071 

0.6396 

-0.2023 

-0.2236 

0 

0 

0 

0 

-0.7071 

-0.6396 

0.2023 

0.2236 

0 

0 

0 

0 

0 

-0.4264 

-0.6068 

-0.6708 

0 

0 

0 

0 

0 

0 

-0.7416 

0.6708 

.4142 

-2.1213 

0 

-0.7071 

0 

0 

0 

0 " 

0 

-2.3452 

-0.4264 

0.6396 

0 

0 

0 

0 

0 

0 

-1.3484 

-2.4271 

0 

0 

0 

0 

0 

0 

0 

1.7889 

0 

0 

0 

0 

0 

0 

0 

0 

-1.4142 

-2.1213 

0 

-0.7071 

0 

0 

0 

0 

0 

-2.3452 

-0.4264 

0.6396 

0 

0 

0 

0 

0 

0 

-1.3484 

-2.4271 

0 

0 

0 

0 

0 

0 

0 

1.7889 , 


Now by using (3. 2), we have 


Pi ' 


f r 3 


6 1.375 + 0.625r ' 

Pi 


4 + r 


0.875 + 0.125r 

<h 


1 + r 


0.125 + 0.625r 



r-4 


-1.375 + 0.125r 


= R'Q t 


= 


~Pi 

r — 2 


-0.2875 + 0.875r 

— 


2r-7 


-1.375 + 0.375/- 

~Pi 





— 


r — 3 


-1.625 + 0.875/- 

~ c h 







U+ij 


v 0.875 + 0.375r y 


\-<h) 


The solution may be written as 


f ~ \ 
Zt 

V^2 j 


f (1.375 + 0.125r, 2.875 -0.875r) + i(0. 125 + 0.625r, 1.625 -0.875r) " 
v (0.875 + 0. 125r, 1.375 -0.375r) + /(-l .375 +0. 125 r, -0.875 - 0.375 r). 


Example 4.2 Conside a simple RLC circuit with fuzzy current and fuzzy sourse in Figure l.The CFSLE for the circuit is as 
follow : 
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(5,, 1,D 
+«( 0 , 1 , 1 ) 



(- 1 . 1 , 1 ) 

+K-2.1.1) 


(10 - 7 . 5 /) z x - (6 - 5 /) Zi = (4 + r, 6 - r) + /(-l + r, 1 - r) 

< 

-(6 - 5 /) Zj + (16 + 3/) z 2 = (-2 + r, - r ) + /(-3 + r, -1 - r) 

Let z, = p x + iq x =(p v p 1 ) + i(q l ,q l ), Z 2 = p^ + iq-, = (p^,q 2 ) + i(q 2 ,q 2 ). The extend matrix is 

^10 0 7.5 0 0 6 0 5 ^ 

0 16 006053 

0 5 10 0 7.5 0 0 6 

530 16 0060 

5 = 

0 6 0 5 10 0 7.5 0 

60530 16 00 

7.5 0 0 6 0 5 10 0 

^ 0 0 6 0 5 3 0 16, 

And we obtain 


^—0.6785 

0.0383 

-0.2174 

0.3857 

-0.0239 

0.3954 

0.4249 0.0629 > 


0 

-0.8876 

0.1659 

0.3010 

0.2183 

-0.0736 

-0.0590 0.1938 


0 

-0.2774 

-0.7365 

-0.1883 

0.0879 

0.2058 

-0.2931 - 

-04573 


-0.3392 

-0.1473 

0.2181 

-0.7970 

0.3458 

0.1555 

0.1656 0.0966 


0 

-0.3382 

0.0622 

-0.2472 

-0.8283 

-0.0528 

0.3188 - 

0.1839 


-0.4071 

0.0230 

-0.1699 

0.0013 

0.1107 

-0.8671 

0.1014 - 

0.1747 


-0.5088 

0.0287 

0.2804 

0.0160 

-0.2873 

0.0628 

-0.7580 - 

0.0169 


v 0 

0 

-4730 

-0.1694 

-0.2090 

-0.1146 

-0.1272 0.8213 y 


14.7394 

-1.0177 

-7.1238 

-9.7019 

0 

-13.1281 7.1238 

-3.3923 

0 

-18.0268 

-2.3715 

-3.7792 

-10.7340 

0.7411 

-7.5305 

-4.1354 

0 

0 

-12.0838 

4.9730 

-6.2713 

-4.0395 

5.4090 

12.5773 

0 

0 

0 

-13.8874 

-2.9250 

1.9066 

-4.9705 

-1.0086 

0 

0 

0 

0 

-7.3583 

-0.4347 

-5.9185 

-2.2809 

0 

0 

0 

0 

0 

-11.5309 0.7967 

1.1572 

0 

0 

0 

0 

0 

0 

-4.4908 

-1.8456 

0 

0 

0 

0 

0 

0 

0 

11.3242 


v 
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Now by using (3. 2), we have 


n. ] 


r 4 + r ^ 


' 0.3 164 + 0.0378r N 

Ei 


—2 + r 


0.0347 +0.0307r 

<h_ 


— 1+ r 


0.0708 + 0.0378r 

c h 

-Pi 

= R~ 'Q T 

-3 + r 
r — 6 

= 

—0.2380 + 0.0307 r 
-0.3920 + 0.0378r 

-P 2 

~<h 

V-?2) 


r 

r-l 

v r + 1 j 


-0.0961 + 0.0307 r 
—0.1464+0.0378r 
v 0.1765 + 0.0307 r , 


The solution may be written as 


( ~ \ 
Zi 

\ Z U 


r (0.3 164 + 0.0378r, 0.3920 - 0.0378r) + /(0.0708 + 0.0378r, 0. 1464 - 0.0378r) " 
v (0.0347 + 0.0307 r, 0.0961 - 0.0307r) + t’(-0.2380 + 0.0307r, -0. 1765 - 0.0307 r). 


V. Conclusion 

In this work we presented a model for solving complex fuzzy linear equation Cz = W where C is a crisp complex matrix 
and W is an arbitrary complex fuzzy vector, respectively. The complex fuzzy linear system is converted to a high order 

linear system SX (r) = F(Vj . We use the QR-decomposition of the coefficient matrix S to obtain fuzzy solution of 

complex fuzzy linear systems. In addition, numerical examples showed that our method is feasible to solve this type of 
complex fuzzy linear systems. 
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